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ABSTRACT 



We present precision radial velocity (RV) data that reveal a multiple exo- 
planet system orbiting the bright nearby G5V star 61 Virginis. Our 4.6 years of 
combined Keck/HIRES and Anglo-Australian Telescope precision RVs indicate 
the hitherto unknown presence of at least three planets orbiting this well-studied 
star. These planets are all on low-eccentricity orbits with periods of 4.2, 38.0, 
and 124.0 days, and projected masses (Msini) of 5.1, 18.2, and 24.0 M ffi , respec- 
tively. Test integrations of systems consistent with the RV data suggest that the 
configuration is dynamically stable. Depending on the effectiveness of tidal dis- 
sipation within the inner planet, the inner two planets may have evolved into an 
eccentricity fixed-point configuration in which the apsidal lines of all three planets 
corotate. This conjecture can be tested with additional observations. We present 
a 16-year time series of photometric observations of 61 Virginis, which comprise 
1194 individual measurements, and indicate that it has excellent photometric 
stability. No significant photometric variations at the periods of the proposed 
planets have been detected. This new system is the first known example of a G- 
type Sun-like star hosting a Super-Earth mass planet. It joins HD 75732 (55 Cnc), 
HD 69830, GJ 581, HD 40307, and GJ 876 as a growing group of exoplanet sys- 
tems that have multiple planets orbiting with periods less than an Earth-year. 
The ubiquity of such systems portends that space-based transit-search missions 
such as KEPLER and COROT will find many multi-transiting systems. 

Subject headings: astrobiology - planetary systems - stars: individual (61 Vir) 
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1. Introduction 

Over 400 extrasolar planets are now known. The majority have been discovered by 
using precision radial velocities (RVs) to detect the reflex barycentric motion of the host 
star. We have had a large sample of over 1000 nearby stars under precision RV survey 
for the past 13 years at Keck and for the past 11 years by the Anglo- Australian Planet 
Search (AAPS) at the Anglo- Australian Telescope (AAT). Of particular interest are the 
nearest brightest stars, as they allow one to achieve the highest precision on stars worthy 
of follow-up with space-based missions such as Hubble Space Telescope (HST), Spitzer, and 
James Webb Space Telescope (JWST). Of equal importance to stellar apparent brightness 
in RV precision is observing cadence. Keplerian signatures are mostly strictly periodic, 
and the signal-to- noise ratio (S/N) of the detection can be enhanced in the presence of a 
random noise background through increased cadence. However, obtaining adequate cadence 
on large telescopes like the Keck and the AAT is quite difficult, as observing time is limited 
for any group. 

Our groups have been combining observations from both the Keck and AAT on select 
stars in the declination overlap region. One of the target stars is 61 Virginis (61 Vir), a 
very nearby G5V star, only 8.5 parsecs away. This star was put on the Keck program in 
2004 December and added to the AAPS target list about four months later. Over the past 
five years, we have accumulated a total of 206 precision RVs that indicate a system of at 
least three planets orbiting this star. In this paper, we present all of these RV data and 
discuss the planetary system that they imply. 
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2. Basic properties of the host star 61 Vir 



61 Vir (=HD 115617, HR 5019, HIP 64924) is a very bright (V=4.74) and very 



well-studied star. It 



pc ( jPerryman et al- 



ias a spectral type of G5V and lies at a distance of only 8.52±0.05 



19971 ). This star has been characterized by a n umber of studies, w ith 



prope rties of interest to 



(I2005f ) and 



SousaetaL 



planet ary system characterization listed by 



Valenti fc Fischer 



(120081 ). Table [T] summarizes recent determinations of the 
fundamental stellar parameters for 61 Vir. Taken together, the properties of 61 Vir indicate 
that it is an old, inactive star ideally suited for preci sion RV planet se arches. Using 



Ca H+K measurements taken between 1994 and 2006, 



Hall et al. 



(120071 ) found that 61 Vir 



is one of 13 targets for which the observed variability is zero within t he uncertainties. In 29 



observations over 7 seasons, they measure a mean log-R^K °f — 4.93. 



Henry et al. 



(1996) 



also measured activity in this st ar, and they found 



with Hall et al. above. Similarly, 



Wittenmyer et al. 



okR' } 



HK 



-4.96, which is consistent 



20061 ) derived a mean \ogR' HK = — 5.03 



from five years (N = 18) of observati ons at McDona ld Observatory. Our measurement of 



\ogR' HK = —4.95 leads to an e stimate (Wright 



due to stellar surface activity. 



Baliunas et al. 



2005) of 1.5 ms 1 for the expected RV jitter 



(119961 ) measured a rot ation period of P = 29 



clays 



or 61 Vir. The age of 61 Vir was es t imate d as 6.3^gj Gyr by 



Valenti fc Fischer 



(l2005f ) and 8.96±££j Gyr by 



Takeda et al. 



(120071 ). In summary, 61 Vir is a nearby, bright, 



solar-type star with physical properties quite similar to our own Sun. 



3. Observations 



The HIR ES spectrometer (|Vogt et al. 



spectrometer (IDiego et al. 



1994J) of the Keck I telescope and the UCLES 



19901 ) of the AAT have been used to monitor 61 Vir. A total of 



126 AAT observations, dating from 2005 April 21 to 2009 August 14, representing a time 
span of 1576 days, have been obtained. The median internal velocity uncertainty for these 
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AAT data is 0.65 ms -1 . A total of 80 Keck observations dating from 2004 December 29 to 
2009 Aug 09, representing a data span of 1685 days, have also been obtained. The median 
internal velocity uncertainty for these Keck data is 0.54 ms _1 . 



Doppler shifts were measured (IButler et al. Ill996l ) by placing an Iodine absorption 



cell just ahead of the spectrometer slit in the converging beam from the telescope. This 
gaseous Iodine absorption cell superimposes a rich forest of Iodine lines on the stellar 
spectrum, providing a wavelength calibration and proxy for the point spread function (PSF) 
of the spectrometer. The Iodine cell is sealed and temperature-controlled to 50 ± 0.1 °C 
(at Keck) and 60 ± 0.1 °C (at the AAT) so that the column density of Iodine remains 
constant. For the Keck planet search program, we operate the HIRES spectrometer at a 
spectral resolving power R « 70, 000 and wavelength range of 3700-8000 A, though only 
the region 5000-6200 A (with Iodine lines) is used in the present Doppler analysis. For the 
AAT program, we typically achieve a spectral resolving power of R ~ 50, 000. Doppler 
shifts from the spectr a are determined with the spectral synthesis technique described by 



O 'Toole et al. 



(120081 ). The Iodine region is divided into ~700 chunks of 2 A each. Each 
chunk produces an independent measure of the wavelength, PSF, and Doppler shift. The 
final measured velocity is the weighted mean of the velocities of the individual chunks. 

Table [2] lists the complete set of 206 RVs for 61 Vir, corrected to the solar system 
barycenter. The table lists the JD of observation center, RV, and internal uncertainty. No 
offset has been applied between the AAT and Keck in this table. The internal uncertainty 
reflects only one term in the overall error budget and results from a host of systematic errors 
from characterizing and determining the PSF, detector imperfections, optical aberrations, 
effects of under-sampling the Iodine lines, etc. Two additional major sources of error 
are photon statistics and stellar jitter. The latter varies widely from star to star, and 
can be mitigated to some degree by selecting magnetically inactive older stars and by 
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time- averaging over the star's unresolved dominant asteroseismological p-mode oscillations. 
Since the single exposures required to reach a required S/N for bright stars like 61 Vir 
are much shorter than the characteristic time scale of low-degree surface p-modes, short 
exposures (ie., <5 minutes) will add an additional noise (or "jitter") component. This 
latter effect was recognized as a potential noise source by the AAPS some time ago, and so, 
since 2005 July, AAPS observations of bright targets like 61 Vir have been extended to th e 



10- 



2008 



5 minu tes periods required to average over these p-mode oscillations (jOToole et al 



20091 ). For most of the past four years, only single exposures of ~7 s at Keck were 
taken of 61 Vir at each epoch, though in 2008 July we began p-mode averaging by combining 
multiple shots of 61 Vir over a 5-10 minutes dwell at each epoch. All observations are then 
further binned on two-hour timescales after being precision Doppler processed. 



4. Photometry 



In addition to our RV observations from Keck I and AAPS, we have acquired 
high-precision photometric observations of 61 Vir during 17 consecutive observing seasons 
from 1993 April to 2009 April with the T4 0.75 m automatic photometric telescope (APT) 
at Fairborn Observatory. Our APTs can detect short-term, low-amplitude brightness 
variability in solar-type stars due to rotational modulation of the visibility of surface 
magnetic activity (spots and plages), as well as longer term variations associated with 
the gro wth and deca y of individual active regions and the occurence of stellar magnetic 



cycles ( jHenry 



19991 ). The photometric observations help to establish whet her observed 



RV va ri ations are cau s ed by stell ar activity or p 



2000ah . 



Queloz et al. 



(1200 lh and 



Paulson et al. 



anetar y reflex motion (e.g., 



Henry et al. 



(120041 ) have presented several examples of 



periodic RV variations in solar-type stars caused by photospheric spots and plages. The 
photometric observations are also useful to search for transits of the planetary companions 
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e.g., iHenrv et al. IbOQObf ) 



The T4 APT is equipped with a precision photometer based on an electromagnetic 
interference (EMI) 9924B bi-alkali photomultiplier tube that measures photon count rates 
successively through Stromgren b and y filters. The APT measures the difference in 
brightness between a program star and nearby constant comparison stars. Our automatic 
tel escopes , photo meters, observing procedures, and data reduction techniques are described 



in 



Henry 



(Il999l ). Further details on the devel opment and operatio n of the automated 



telescopes can be found in iHenryl (11995a 



bj) and 



Eaton et al. 



(|2003j). 



For 61 Vir, we used the two comparison stars HD 113415 (CI, V = 5.58, B — V = 0.56, 
F7 V) and HD 114946 (C2, V = 5.33, B -V = 0.87, G8 III-IV). The individual Strdmgren 
b and y differential magnitudes have been corrected for differential extinction with 
nightly extinction coefficients and transformed to the Stromgren system with yearly mean 
transformation coefficients. Since 61 Vir lies at a declination of —18°, the photometric 
observations from Fairborn were made at air mass 1.6-1.8. Therefore, to maximize the 
precision of the measurements, we combined the Stromgren b and y differential magnitudes 
into a single (b + y)/2 passband. We also computed the differential magnitudes of 61 Vir 
with respect to the mean brightness of the two comparison stars. Because we are interested 
only in variability timescales of days to weeks, we have normalized the final differential 
magnitudes so that the mean brightness of each observing season is equal to zero. This 
effectively removes any long-term brightness variability in the two comparison stars as well 
as in 61 Vir. (However, the standard deviation of the seasonal mean differential magnitudes 
was only 0.00033 mag before normalization.) 

A total of 1194 normalized differential magnitudes from 17 observing seasons are 
plotted in the top panel of Figure [U The data scatter about their mean with a standard 
deviation a = 0.00196 mag, which provides an upper limit to possible brightness variation 
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in 61 Vir. A Lomb-Scargle periodogram of the photometric measurements is shown in 
the bottom panel of Figure [T] and reveals no significant periodicities within the data. We 
computed least-squares sine fits for the three RV planet candidates described below; the 
semi-amplitudes of the lightcurve fits were 0.00016 ± 0.00007, 0.00011 ± 0.00008, and 
0.00014 ± 0.00007 mag for the three periods 4.215, 38.012, and 123.98 days, respectively. 
We conclude that the photometric constancy of 61 Vir supports planetary reflex motion as 
the cause of the RV variations described in the next section. 



5. The Planetary System Orbiting 61 Vir 

The combined RV data from the AAPS and Keck telescopes show a root-mean-square 
(rms) scatter of 4.1ms _1 about the mean velocity. This includes an offset of 0.895 
ms _1 between the two telescopes (Keck - AAT) that was left as a free parameter and 
emerged from the three-planet Keplerian fit. This rms significantly exceeds both the 
scatter due to the underlying precisions of both the Keck and AAPS Doppler measurement 
pipelines and the scatter expected in this star due to its predicted 1.5 ms" 1 level of stellar 
jitter. Figure [2] shows the combined RV data set. 

Figure [3] (top panel) shows the circular periodogram of the combined RV data set. 
This is the power at any period associated with fitting circular orbits to the data. Power 
at each sampled period is proportional to the relative improvement in the fit quality for a 
circular model versus a constant velocity model — that is, the relative drop in xl- Figure 
[3] (bottom pane l) shows the spectral window or power spectrum due to the sampling times 



( iDeeming 



19751 ). This spectral window illustrates spurious power that can appear in the 



data merely from the sampling times alone. 

The periodogram of the RV data (top panel of Figure [3]) shows a number of significant 
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Fig. 1. — Differential photometry of 61 Vir (top panel). Periodogram of the photometry 
(bottom panel). 
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Fig. 2. — Relative radial velocities of 61 Vir. Velocities obtained at the AAT are shown in 
red, and those from Keck are shown in blue. 
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signals, with the strongest peak occurring at a period of P = 38.13 days. The false 
alarm probability (FAP) of this peak is estimated (adopting the procedure described in 



Cumming 



(|2004Ji ) to be FAP~ 3.8 x 1(T 16 . The horizontal lines in Figure [3] (top panel) 
and all similar figures below indicate (top to bottom) the 0.1%, 1%, and 10% FAP levels. 
This P = 38.13-day signal, furthermore, lies far from the periods favored by the sampling 
window (Figure [31 bottom panel), which produces spurious power at periods near 357, 191, 
95.5, 652, 29.6, and 110 days. (Signals near these periods would naturally attract suspicion 
as being artifacts of the lunar and yearly periods on telescope scheduling.) 

The mass of t he ho st star is assumed to be 0.95 M , the isochrone mass of 



Valenti fc Fischer 



(120051 ). Based on the periodogram, we then fit a planet of mass 
Msinz = 15M e and P = 38.20 days on a circular orbit to the RV data. The presence 
of this planet (with RV semi-amplitude K = 2.98 ms" 1 ) reduces the rms scatter of the 
velocity residuals to 3.39 ms -1 . Figure H] (top panel) shows the periodogram of the 
residuals to the one-planet fit which has strong peaks at P = 124 days and P = 4.21 days. 
The P = 124-day signal has FAP~ 6.3 x 10~ 12 and can be modeled with a companion with 
K = 3.36 ms -1 and Msini = 24M e . The addition of this planet further reduces the rms 
scatter to 2.78 ms -1 . 

Figure H] (bottom panel) shows the periodogram of the residuals to the two-planet fit 
which has a highly significant peak at a period of P = 4.215 days. As a consequence of the 
fact that the outer two planets, once removed, have significantly reduced the variance in the 
RV data, the FAP of the 4.215-day signal in the residuals is very small, FAP~ 1.1 x 10~ 22 . 
We are thus confident that this periodicity is real, and we ascribe it to the presence of a 
P = 4.2149 days, Msini = 5.1M® companion. There is another peak at 1.3 days that arises 
as an alias of the 4.2-day periodicity beating with the nightly sampling. The three-planet 
model has an rms scatter of 2.17 ms -1 . 
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Fig. 3. — Top panel: circular periodogram of the combined RV data set for 61 Vir. Bottom 
panel: power spectral window of the combined RV data for 61 Vir. 
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Fig. 4. — Top panel: circular periodogram of the one-planet residuals of the combined RV 
data set for 61 Vir. Bottom panel: circular periodogram of the two-planet residuals of the 
combined RV data set for 61 Vir. 
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Given the three-planet model, we can look either for solutions in which the planetary 
orbits are circular or solutions where the eccentricities are allowed to float. Inclusion of 
eccentricities provides only a modest improvement to the orbital fit, and we conclude that 
a significant amount of additional Doppler velocity monitoring will be required to improve 
measurement of the eccentricities. For a detailed discussion o f how additional measurements 
can reduce uncertainties in orbital elements, see iFordl (120051 ). 



In Tables |3] and HI we present our best-fit versions of the system under the assumption 
of circular orbits (Table [3D and with the additional degrees of freedom provided by fully 
Keplerian trajectories (Table HJ). For the orbital fits, we assume i = 90° and Q = 0°. 
We have verified that the inclusion of planet-planet gravitational interactions in the fit 
are unnecessary. Most of our modeling involves simply adding Keplerians. However, for 
cases where we expect significant gravitational interactions between companions, we carry 
out a more detailed modeling that involves using Hermite fourth order approximations to 
the planets' trajectories, accurately characterizing such gravitational interactions. If the 
Hermite fourth order calculations do not differ significantly from simple summed Keplerians, 
as is the case here, we conclude that planet-planet interactions are not necessary. 

Uncertainties ar e based on 1 00 bo otstrap trials for which we follow the procedure 



in Section 15.6 from 



Press et al. 



(119921 ). We take the standard deviations of the fitted 
parameters to the bootstrapped RVs as the uncertainties. The fitted mean anomalies are 
reported at epoch JD 2453369.166. Our fitting was carried out with the publicly available 



Systemic Console (IMeschiari et al 



20091 ) 



Figure [5] (top panel) shows the power spectrum of the velocity residuals for the 
three-planet e = fit. There is a significant peak near 94 days with FAP~ 0.0003. Although 
using this period in a four-planet circular fit results in a significant improvement in xt, 
the rms decrease from 2.17 ms 4 for the three-planet model to 2.00 ms" 1 is not significant. 
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Additionally, the fitted K = 1.44 ms -1 is significantly smaller than the scatter around the 
model. But perhaps most significantly, this fourth peak almost exactly corresponds to a 
peak in the window function of our data (Figure [31 lower panel), making any association of 
it with a real planet questionable. 

The chi-squared of our three-planet circular fit is 13.03 and results in a fit with an rms 
of 2.17 ms -1 and estimated stellar jitter of 2.06 ms -1 . That estimate of the stellar jitter is 
the jitter required to bring the chi-squared of the fit down to 1.0. Thus, if the stellar jitter 
is 2.06 ms -1 , our three-planet fit is essentially perfect. The difference between a true stellar 
jitter of 2.06 ms -1 and our estimated value of 1.5 ms -1 is negligible, given the accuracy of 
such jitter estimates. Moreover, there are likely to be other planets in this system such as 
the 94-day signal in the top panel of Figure [5j Adding more planets to the model would 
further reduce the stellar jitter component of the fit. But the present data set does not 
support adding another Keplerian to the model. 

Figure [5] (bottom panel) also shows the power spectrum of the velocity residuals for the 
eccentric three-planet fit. The tallest peak is near 10.9 days, but it is not significant with 
a FAP near 20%. Also, note that the rms for the eccentric three-planet fit of 2.09 ms _1 is 
not significantly different from the 2.17 ms -1 for the circular three-planet fit. At this point, 
all of these concerns indicate that many more observations are needed both to constrain 
the eccentricities of the first three planets and to pursue the prospect of a potential fourth 
planet in the system. 

To summarize up to this point, we have shown that the combined Keck and AAT RV 
data show strong evidence for a three-planet system in orbit about the G5V star 61 Vir. 
Figure [6] shows the barycentric RV of the star computed with the eccentric three-planet 
model. Figure [7] shows the barycentric RV of the star computed with the eccentric model 
due to each individual companion in the system. In each panel, the velocities are folded 
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Fig. 5. — Top panel: circular periodogram of the three-planet residuals, assuming circular 
orbits, of the combined RV data set for 61 Vir. Bottom panel: circular periodogram of the 
three-planet residuals, with fitted eccentricities. 
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at the period of each corresponding planet. Points marked as open symbols denote 
observations that are suspect due to poor observing conditions that led to low S/N in the 
spectra. All such points were previously noted in the observing logs before any analysis 
occurred. Regardless of such suspicions, all observations were used in all the analyses. 



6. Dynamical Analysis 



For any multiple-planet system, long-term orbital stability is a matter of both interest 
and concern. Even precisely characterized, seemingly stable configurations such as the Solar 
System are subject to chaotic planet-planet interactions that can lead to orbit-crossing, 



collisions , and ejections o n time scales shorter than t 



see, e.g., 



Laskar 



fl2008h : 



Batygin &: Laughlin 



re stellar lifetime (for recent t reatments 



2008|); 



Laskar fc Gastineau 



(120091 ) I Given 



the uncertainties that adhere to the orbital elements, the effort of generating a dynamical 
characterization of a system such as 61 Vir b-c-d, is somewhat akin to producing detailed 
maps of sand dunes. Nonetheless, it is useful to verify, via numerical integration, whether 
the planetary configurations listed in Tables [3] and H] are dynamically stable. Such an 
analysis is particularly useful in giving rough bounds on the allowed coplanar inclinations 
relative to the line of sight to Earth. For simulations that do not include the effects of 
tidal dissipation, we use fitted Newtonian parameters for the initial states for long-term 
i ntegrations . Wit h a time step of 0.1 day, we use the MERCURY integration package 



( 1 Chambers 



19991 1 for th e simulations. We i n clude the first order post-Newtonian term in 



the star's potential as in 



Lissauer fc Rivera 



(120011 ). 



If we assume the three orbits are initially circular, with periods, masses, and mean 
anomalies given in Table [3J we find that the system is stable for at least 365 Myr. 
Additionally, assuming the system to be coplanar, if we set the inclination to the sky plane 
to various values from % = 90° all the way down to i — 1° and perform a Newtonian fit for 
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Time (JD - 2450000.0) 

Fig. 6. — Top panel: observed RVs of 61 Vir (red points for AAT observations and blue 
for those from Keck) and model RV of 61 Vir due to the eccentric three-planet fit listed in 
Table H] (black curve). The fitted relative offset of 0.895 ms" 1 has been applied between the 
two data sets from the two telescopes. Bottom panel: residual velocities remaining after 
subtracting the eccentric three-planet model from the observations (again, red is for AAT, 
and blue is for Keck). 
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Fig. 7. — Top panel: radial velocity of 61 Vir due to planet b folded at 4.215 days. Center: 
radial velocity of 61 Vir due to planet c folded at 38.021 days. Bottom panel: radial velocity 
of 61 Vir due to planet d folded at 123.01 days. In each panel, the effect of the other two 
planets has been subtracted out. The curves represent the model velocities due to each 
respective planet. The AAT observations are shown in red, and the Keck velocities are 
shown in blue. Open symbols represent observations that may be suspect due to low S/N 
based on observing log notes. 
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the other 11 parameters (three parameters per planet plus the two velocity offsets), xt does 
not change significantly from the nominal % = 90° fit. We also find the % = 1° fit to result 
in a system that is stable for at least 50 Myr. For this inclination, the fitted masses exceed 
1.0, 3.3, and 4.5 Mj up . This system is stable because of the small eccentricities. Thus, 
under the assumption that the system is coplanar and the orbits are (nearly) circular, we 
cannot place a lower bound on the inclination of the system. 

The parameters of the floating-eccentricity version of the 61 Vir system given in Table H] 
were also used as the initial input conditions for a 10 7 yr simulation. The three-planet 
configuration remained stable for the full duration of this simulation. Figure [H] shows the 
evolution of the planets' semi-major axes and eccentricities over the first 250,000 years of 
the simulation. The planets interact with each other on secular timescales, and the inner 
planet experiences a large excursion in eccentricity (up to e ~0.5) but is never disrupted. 

With its period of 4.215 days, the innermost planet in the system is likely subject to 
significant tidal dissipation in its interior. Tidal damping, in fact, may have played an 
important role in the evolutionary history of this multi-planet system. An investigation 
of this process requires a dynamical theory that couples dissipative orbital evolution with 
planet-planet gravitational perturbations among the individual system members. While the 
coupling between dissipative and secular processes can be complex, there exists an avenue 
toward insight regarding the overall dynamical state of the system. In particular, it may be 
possible that the dissipative properties of the innermost planet can be inferred. 

The key idea is that, under the influence of tidal dissipation, multiple-planet systems 
approach stationary configurati ons. These so-called tida l fixed points are discus sed in the 



contex t of extrasolar planets by 



Wu k Goldreich 



(120021 ) and more generally by 



Mardling 



( 120071 ). Fixed-point configurations are characterized by either parallel or anti-parallel 



alignment of the apsides, and simultaneous precession of the apsidal lines. Additionally, the 
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Fig. 8. — Top panel: semi-major axes vs. time for the simulation using the parameters from 
Table H] as starting conditions. Bottom panel: eccentricities vs. time for the same simulation. 
For both panels, planet b is in red, c is in blue, and d is in green. 
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eccentricities of such systems are well determined and are not subject to significant time 
variations. 

Participation by the system in a fixed-point configuration hinges on the tidal quality 
factor, Q, of planet b. Assuming that the system did not form in a stationary configuration, 
and that the dissipation of planetary orbital energy has been dominated during the past 
by tides raised on the innermost plane t , it ca n be shown that it takes ~ 3 circularization 



timescales, r, given by 



Yoder fc Peale 



(1981) 



Q 



m 



X 



21n \ k j \ MJ \R, 

where k is the plan etary Love numb er and R is the inner planet's radius, for the system to 



become stationary (IMardling 



20071 ) . Therefore, assuming R = 2.45i? e and k = 0.34, we 



expect that it would take ~ 700, 000 x Q years for the system to arrive at the fixed point. 
We have co nfirmed this, by perform ing a numerical integration , using the tidal friction 



formalism of 



Eggleton et al. 



dlSaSP and 



Mardling fc Lin 



(120021 ) . and the eccentric fit listed 



in Table H] as initial conditions. Indeed, libration of planet b around the fixed point begins 
after ~ 850, 000 x Q years, in rough agreement with theory. It is further interesting to 
note that the best-fit eccentricities for planets c and d are in fact already very close to the 
fixed-point values, and their perihelia are already set in a librating mode at the beginning 
of the integration. 

Thus, given the star's multi-billion year age, we can make the following claim: if planet 
b has a Neptune-like Q (~ 30,000), then the system is probably not stationary. However, 
if the planet has a characteristic terrestrial Q (~ 100), then we can expect the system to 
be at a fixed point. If the system is at a fixed point, and the eccentricities are modest, 
we can use modified Laplace-Lagrange secular theory to relate the planets' eccent r icities 



to each-other. The three-planet fixed point equations presented in 



Batygin et al. 



(120091) 



yield the relationships: e c ~ 5.11e&, and ~ 1 .TJe h . Thus, by a precise determination of 



-23- 



fixed-point eccentricities, we can not only determine the dynamical state of the system, but 
also begin to characterize the nature of planet b. 



7. Discussion 



We have detected three low-mass planets orbiting the nearby star 61 Vir. With a RV 
semi-amplitude K = 2.09 ±0.23 ms" 1 , the Msini = 5.10 M® inner planet 61 Virb is among 
the very lowest amplitude companions yet detected using the Doppler velocity technique. 
To date, the only annou nced planets with sm aller amplitudes have been Gliese 581e, with 



K 



.85 ± 0.23 ms" 1 fMavor et al. 



ms 1 f Mavor et al 



ms 1 (ILovis et al. 



2009bJ), and HD 40307b, with K = 1.97 ± 0.11 



2009aJ). The inner planet in the HD 69830 system, with K = 2.20 ±0.20 



20061 ) presents a slightly larger K. Interestingly, all four of the lowest 



amplitude planets known are members of multiple-planet systems containing three or more 
Neptune or super-Earth class objects on orbits that would be ascribed to the inner zone of 
the terrestrial region were they in our own solar system. 

Given the observed multiplicity, negligible jitter, and extraordinary photometric 
stability associated with the 61 Vir planetary system, it is natural to speculate on the 
existence and detectability of additional planets. As 61 Vir is a naked-eye visible solar-type 
star observable from both hemispheres, it has been the target of many plan et-search effort s 
i n the past 25 years. One of the first precision RV searches for planets, by 



Walker et al. 



( 119951 ). observed 61 Vir for approximately 12 years. They present a figure indicating 



rough upper limi t s of 1 -2 Jupiter masses in the period range between one and ten years. 



Cumming et al. 



(Il999l ) included 61 Vir in their determination of companion limits from 



11 yr of the Lick planet search. They computed a 99% upper limit on the velocity amplitude 
K — 19 ms -1 and n oted a possible candidate companion with a period of 9.8 yr and 



K ~ 15 ms 1 . Later, 



Wittenmyer et al. 



(120061 ) combined 23 years of velocities for 61 Vir, 



-24- 



including those data from IWalker et al. I (119951 ) and from the McDonald Observatory 
long-term planet search program, to derive companion-mass upper limits of 1.16 Mj up at 
3 AU and 1.69 Mj up at 5.2 AU. From these studies, and the lack of any significant trends in 
the RV data, it would appear that there are no planets in the 61 Vir system with masses 
substantially greater than that of Jupiter. Undetected sub- Jovian mass planets may be 
present in orbits with periods longer than the extant high-precision velocity data. 

A further intere sting clue to the nature of the 61 Vir system was recently provided by 



Tanner et al. 



(120091 ). who detected an infrared excess at 160/z using Spitzer observations. 
The excess 160/i flux implies the presence of a dust disk that is continuously replenished 
by collisions within a cold Kuiper-Belt-like disk of planetesimals surrounding the star. 



Tanner et al. 



(120091 ) also suggest that the disk is resolved at 70/x (G. Bryden et al. 2009, 
in preparation). Assuming that the emitting grains are black bodies, the disk spans an 
annulus from Ri nner = 96±5 AU to -Router = 195±10 AU, with the dust temperature at 
55-45 K (inner to outer). An improved model of the disk, however, can be achieved by 
assuming the emission arises from silicate grains with size 0.25//. In this case, the disk 
spans .R irm er = 120±20 AU to .R ou ter = 220±10 AU with a dust tempe rature of 24-19 K. An 



excess at 70 /x was also detected with Spitzer by 



Trilling et al. 



(120081 ) . 



The 61 Vir syst em is consistent with the population of planets postulated by 



Mayor et al. 



(j2009al ). who inferred that fully one-third of solar-type stars in the immediate 
galactic neighborhood are accompanied by Neptune- (or lower) mass planets having orbital 
periods of 50 days or less. This putative population was largely unexpected, and has 
spurred a good deal of recent theoretical work geared to explain its existence. We note that 
the habitable zone of 61 Vir is located in the vicinity of a 300-day orbit. At this separation 
from the parent star, a 2M e planet induces a RV semi-amplitude K = 0.2 ms _1 . Assuming 
a combined stellar jitter and median internal uncertainty of 1.5 ms" 1 , this is a factor of 
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about 30 greater than the 0.05 ms -1 precision required to make a 4cr measurement of a 
such a small signal. If precision scales as the square root of the total number of Doppler 
measurements, a total of 30 2 or 900 velocities (an additional 700 more) would allow for the 
detection of such a planet with a S/N of 4. 

With a period of only 4.215 days (a =0.05 AU), 61 Virb may transit its star. The a 
priori transit probability is ~ 9% (assuming — 1 Rq). Due to the low mass of 61 Virb, a 
successful transit observation would allow a detailed physical probe of a fundamentally new 
class of planet. Currently, the low est mass transiting planet is CoRoT-7b, with a radius of 



1.68 ± 0.09 .R© flLeger et al. 



20091 ). CoRoT-7b orbits a G9V star and has a transit depth of 
only AF/F = 3.4 x 10~ 4 . Its mass is estimated by Queloz et al. (2009) to be 4.8 ± 0.8M e . 
By comparison, if 61 Vir b is transiting, sin i ~ 1, and hence it would have a mass of 
5.10 M®. Assuming that the planet migrated inward from beyond th e ice line of 61 Vir's 



Fortney et al. 



(|2007|) 



protoplanetary disk, it is likely made mostly of water. The models of 
indicate a 2.45 Rq radius for such a planet, leading to a central transit depth of ~ 0.54 
millimagnitudes. Such a transit could be detected from space using, for example, the Warm 
Spitzer platform. 



Of the currently known planets discovered by the RV technique with Msini < 10M f 



OOj) 



anets; 6 host stars), half orbit M dwarf stars. Those are GJ 58 



2009bf ). GJ 876d flRivera et al. 



2005|), and GJ 176b ( 



those planets orbit K stars: HP 7924b ( 



20091 ). and HD 40307b,c,d flMavor et al. 



oward et al. 



brveille et al. 



c,d,e (Mayor et al. 



20091). The rest of 



2009|), HD 181433b flBouchv et al. 



2009al ). This is of course a selection effect, since 



late-type stars have lower masses and hence a planet of a given mass would produce a 
larger RV signal. 61 Vir, therefore, is a unique new system in that not only does it host 
multiple low-mass planets, but also it is the first G-type Sun-like star found through the 
RV technique to host a planet with Msinz < 10M e . 



-26- 



The detection of the three low- if planets reported here around this previously 
well-studied, bright, nearby star was made possible by the combined cadence and high 
precision of the AAPS and Keck surveys. As cadence and time bases grow in the RV 
monitoring of chromospherically quiet, nearby stars, complex planetary systems, like 
55 Cnc, and now 61 Vir, are becoming increasingly common. This is only our first 
reconnaissance of this fascinating and quite nearby system. As more RV data are collected, 
the orbital ephemerides of the planets will become better determined, and more planets will 
probably be revealed. 

The 61 Vir system joins a growing class of exoplanet systems that have multiple planets 
orbiting with periods less than an Earth-year. Other examples are HD 75732 (55 Cnc), 
HD 69830, GJ 581, HD 40307, and GJ 876. The increasing frequency of such systems 
portends that space-based transit surveys such as COROT and KEPLER will find many 
multi-transiting systems. 
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Table 1. Stellar Parameters for 61 Vir 



Parameter 


Value 


Reference 




Spec. Type 


G5V 


Cenarro ct al. 


(2007) 


iviass yivi ) 


U ' yt, -0.03 


Valenti & Fischer 


(2005) 


Radius (Rq) 


n Q49+ 034 

u.»^_ 02g 
0.963±0.011 


Takcda ct al. 
Valenti & Fischer 


(2007) 
(2005) 




0.98±0.03 


Takcda ct al. 


(2007) 


Luminosity (Lq) 


0.805±0.028 


Valenti & Fischer 


(2005) 




n cri/i _i_n nn^ 


Sousa ct al. 


(2008) 


Distance (pc) 


8.52±0.05 


Pcrrvman et al. 


(1997) 


V sin % (km s ^ ) 


2 2 


Valenti & Fischer 


(2005) 






Dcsidera ct al. 


(2006) 


logi^ K 


-4.93 


Hall ct al. 


(2007) 






Henry ct al. 


(1996) 




-5.03 


Wittcnmycr ct al. 


(2006) 




-4.95 


This work 


P ro t (days) 


29 


Baliunas ct al. 


(1996) 


Age (Gyr) 


fi q+3.3 
°-°-3.1 


Valenti & Fischer 


(2005) 


[Fe/H] 


B - 9b -3.08 
0.05 


Takcda ct al. 
Valenti & Fischer 


(2007) 
(2005) 




-0.01 


Cenarro ct al. 


(2007) 


Teff (K) 


5571 


Valenti & Fischer 


(2005) 




5531 


Cenarro ct al. 


(2007) 




5577±33 


Ecuvillon ct al. 


(2006) 




4.47 


Valenti & Fischer 


(2005) 




4.31 


Cenarro ct al. 


(2007) 




4 45+0.04 
4.34±0.03 


Takcda ct al. 
Ecuvillon ct al. 


(2007) 
(2006) 
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Table 2. Radial Velocities for 61 Vir 



JD 


RV 


Error 


Observatory 


(-2450000) 


(ms 4 ) 


(ms 4 ) 




3482.96440 


2.10 


0.73 


A 


3483.96117 


5.56 


0.60 


A 


3485.07996 


3.54 


0.59 


A 


3486.04514 


0.25 


0.48 


A 


3488.13906 


2.37 


0.59 


A 


3489.08696 


3.00 


0.61 


A 


3516.04043 


3.76 


0.71 


A 


3521.94298 


1.94 


0.96 


A 


3569.90141 


-4.25 


0.83 


A 


3570.94617 


-3.47 


0.57 


A 


3571.94048 


-1.64 


0.63 


A 


3573.87030 


-2.32 


0.64 


A 


3575.87530 


1.93 


0.61 


A 


3576 89023 


4.07 


0.55 


A 


3577.86776 


3.57 


0.67 


A 


3578.91377 


-1.40 


0.57 


A 


3840.21603 


4.18 


0.77 


A 


3841.16188 


7.05 


0.65 


A 


3843.13891 


10.21 


0.62 


A 


3844.08319 


6.07 


0.66 


A 
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Table 2 — Continued 



JD 


RV 


Error 


Observatory 


(-2450000) 


(ms 4 ) 


(ms 4 ) 




3937.93317 


-2.73 


0.72 


A 


3938.93436 


0.21 


0.64 


A 


3943.88525 


-2.85 


0.53 


A 


3944.88569 


-5.13 


0.47 


A 


3945.88638 


-6.74 


0.55 


A 


3946.88960 


-3.68 


0.55 


A 


4111.21051 


7.14 


0.77 


A 


4112.21939 


6.16 


0.60 


A 


4113.22990 


4.11 


0.75 


A 


4114.26145 


3.12 


0.90 


A 


4119.24597 


3.94 


0.74 


A 


4120.22002 


3.31 


0.70 


A 


4121.21210 


1.44 


0.56 


A 


4123 23541 


-3.13 


0.63 


A 


4126.19788 


-3.12 


0.65 


A 


4127.20219 


-3.88 


0.62 


A 


4128.20519 


-2.82 


0.73 


A 


4129.20022 


0.69 


0.52 


A 


4130.19130 


-4.28 


0.53 


A 


4131.20082 


-5.25 


0.53 


A 
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Table 2 — Continued 



JD 


RV 


Error 


Observatory 


(-2450000) 


(ms 4 ) 


(ms 4 ) 




4132.20696 


-4.95 


0.73 


A 


4135.19894 


-7.04 


0.67 


A 


4136.21567 


-4.83 


0.56 


A 


4137.21358 


-0.90 


0.49 


A 


4138.19311 


-1.62 


0.82 


A 


4139.18689 


-3.61 


0.72 


A 


4141.21764 


-0.05 


0.71 


A 


4142.20644 


0.61 


0.51 


A 


4144.14374 


-0.35 


0.57 


A 


4145.17237 


2.53 


0.58 


A 


4146.19289 


5.15 


0.55 


A 


4147.21096 


-0.63 


0.55 


A 


4148.24103 


-3.18 


0.57 


A 


4149 18109 


2.16 


0.63 


A 


4150.21285 


3.30 


0.62 


A 


4151.22296 


0.56 


0.64 


A 


4152.09352 


0.80 


0.73 


A 


4152.25396 


-2.16 


0.54 


A 


4153.17187 


2.83 


0.68 


A 


4154.09474 


4.38 


0.57 


A 
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Table 2 — Continued 



JD 


RV 


Error 


Observatory 


(-2450000) 


(ms _1 ) 


(ms 4 ) 




4154.27050 


3.11 


0.66 


A 


4155.06551 


2.80 


0.59 


A 


4155.27322 


2.43 


0.64 


A 


4156.18613 


0.44 


0.42 


A 


4157.16889 


2.37 


0.89 


A 


4158.18706 


3.65 


0.88 


A 


4223.12188 


8.13 


0.54 


A 


4224.15218 


6.77 


0.72 


A 


4225.08458 


10.59 


0.62 


A 


4226.02062 


14.20 


0.85 


A 


4227.02043 


8.15 


0.68 


A 


4252.97759 


-0.11 


0.80 


A 


4254.01661 


2.11 


0.70 


A 


4254 91868 


5.69 


0.85 


A 


4255.98460 


6.39 


0.81 


A 


4257.07180 


1.58 


0.83 


A 


4333.86604 


-0.54 


1.34 


A 


4334.86018 


3.40 


0.85 


A 


4335.85304 


8.92 


0.80 


A 


4336.84663 


7.03 


0.71 


A 
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Table 2 — Continued 



JD 


RV 


Error 


Observatory 


(-2450000) 


(ms _1 ) 


(ms 4 ) 




4543.06967 


0.10 


0.71 


A 


4550.11293 


-2.14 


0.85 


A 


4551.09816 


-2.56 


0.86 


A 


4552.14240 


-5.26 


0.93 


A 


4553.10980 


-2.94 


1.00 


A 


4897.20941 


-4.86 


1.06 


A 


4900.19861 


-0.16 


1.03 


A 


4901.16839 


0.60 


0.80 


A 


4902.21654 


-4.11 


0.99 


A 


4904.20346 


-1.41 


1.53 


A 


4905.27335 


3.23 


0.80 


A 


4906.22330 


-0.55 


0.93 


A 


4907.21266 


0.28 


0.85 


A 


4908 21531 


7.41 


0.75 


A 


5013.83842 


-0.28 


0.61 


A 


5015.83950 


-1.99 


0.68 


A 


5017.84841 


3.13 


0.64 


A 


5018.90225 


0.90 


0.70 


A 


5018.95764 


4.79 


1.09 


A 


5019.99765 


-0.78 


0.73 


A 
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Table 2 — Continued 



JD 


RV 


Error 


Observatory 


(-2450000) 


(ms _1 ) 


(ms 4 ) 




5020.83414 


0.35 


0.60 


A 


5021.88251 


3.15 


0.79 


A 


5022.89587 


3.61 


0.74 


A 


5023.87589 


-0.93 


0.89 


A 


5029.86175 


-0.04 


0.59 


A 


5030.97802 


1.93 


0.50 


A 


5031.84137 


-0.95 


0.60 


A 


5032.92564 


-3.75 


0.53 


A 


5036.86026 


-7.10 


0.63 


A 


5037.84849 


-6.51 


0.71 


A 


5040.85306 


-7.42 


1.04 


A 


5041.95419 


-8.68 


0.84 


A 


5043.84520 


-3.45 


0.52 


A 


5044 84224 


-6.79 


0.49 


A 


5045.83918 


-7.89 


0.62 


A 


5046.91766 


-7.33 


0.64 


A 


5047.86724 


-1.08 


0.47 


A 


5048.85263 


-3.21 


0.47 


A 


5049.86609 


-7.39 


0.46 


A 


5050.87397 


-3.02 


0.68 


A 



-39- 



Table 2 — Continued 



JD 


RV 


Error 


Observatory 


(-2450000) 


(ms 4 ) 


(ms 4 ) 




5051.86595 


-2.08 


0.53 


A 


5052.88848 


-2.74 


0.66 


A 


5053.88182 


-3.97 


0.54 


A 


5054.86448 


-2.46 


0.56 


A 


5055.85374 


2.20 


0.53 


A 


5058.89140 


-1.51 


0.59 


A 


3369.16663 


1.69 


0.73 


K 


3425.07667 


7.30 


0.61 


K 


3479.96032 


0.46 


0.58 


K 


3546.78158 


2.94 


0.44 


K 


3547.76796 


3.03 


0.36 


K 


3548.80045 


-2.27 


0.37 


K 


3549.80104 


-0.03 


0.37 


K 


3550 82365 


-0.08 


0.36 


K 


3551.80855 


-1.31 


0.39 


K 


3576.75813 


6.56 


0.44 


K 


3602.72508 


5.09 


0.49 


K 


3754.09220 


0.80 


0.59 


K 


3776.17013 


4.08 


0.58 


K 


3777.08603 


1.79 


0.56 


K 



-40- 



Table 2 — Continued 



JD 


RV 


Error 


Observatory 


(-2450000) 


(ms _1 ) 


(ms 4 ) 




3806.94545 


-1.23 


0.58 


K 


3926.74769 


5.03 


0.47 


K 


3927.78252 


0.49 


0.42 


K 


3931.82248 


-1.03 


0.44 


K 


4131.14294 


-7.91 


0.65 


K 


4246.81083 


0.59 


0.73 


K 


4247.94552 


-2.04 


0.72 


K 


4248.81966 


-4.38 


1.05 


K 


4250.80828 


0.78 


1.09 


K 


4251.81149 


1.78 


1.02 


K 


4255.75898 


7.13 


0.77 


K 


4277.73265 


-2.06 


0.51 


K 


4278.73377 


-5.19 


0.86 


K 


4279 73075 


-0.88 


0.83 


K 


4280.74026 


-0.24 


1.00 


K 


4285.76472 


-6.15 


0.63 


K 


4294.73843 


-5.39 


0.51 


K 


4304.73624 


4.05 


0.60 


K 


4305.73658 


4.59 


0.61 


K 


4306.80060 


3.51 


0.59 


K 



-41 - 



Table 2 — Continued 



JD 


RV 


Error 


Observatory 


(-2450000) 


(ms _1 ) 


(ms 4 ) 




4307.73485 


-2.50 


0.48 


K 


4308.73551 


2.89 


0.56 


K 


4309.73851 


2.20 


0.77 


K 


4310.73545 


-1.05 


0.84 


K 


4311.73296 


2.89 


0.84 


K 


4312.73423 


-0.92 


0.53 


K 


4313.73496 


1.04 


0.85 


K 


4314.73608 


-2.67 


0.87 


K 


4318.77353 


1.92 


0.51 


K 


4455.16774 


-0.98 


1.19 


K 


4456.15784 


4.78 


1.06 


K 


4465.13890 


0.74 


1.13 


K 


4491.10891 


7.68 


1.08 


K 


4545 03961 


-3.86 


1.10 


K 


4546.03027 


4.83 


1.21 


K 


4547.03936 


0.14 


1.09 


K 


4600.92698 


8.84 


1.07 


K 


4601.87647 


7.50 


0.97 


K 


4634.78181 


8.36 


0.94 


K 


4666.78904 


-4.45 


0.30 


K 



-42- 



Table 2 — Continued 



JD 


RV 


Error 


Observatory 


(-2450000) 


(ms _1 ) 


(ms 4 ) 




4667.78945 


-1.98 


0.29 


K 


4671.80264 


-2.79 


0.30 


K 


4672.78558 


-1.72 


0.30 


K 


4673.79754 


-2.59 


0.29 


K 


4675.79484 


-2.07 


0.51 


K 


4701.73724 


2.44 


0.36 


K 


4702.73348 


4.02 


0.35 


K 


4703.73110 


-1.75 


0.46 


K 


4819.16638 


1.75 


0.37 


K 


4820.17627 


3.19 


0.47 


K 


4821.14761 


0.74 


0.85 


K 


4822.14603 


-2.32 


0.64 


K 


4823.16898 


0.83 


0.38 


K 


4873 10516 


3.59 


0.48 


K 


4903.14262 


-3.53 


0.45 


K 


4967.98451 


0.49 


0.40 


K 


4968.97445 


0.63 


0.74 


K 


5021.83208 


4.26 


0.32 


K 


5022.84171 


5.41 


0.30 


K 


5023.74658 


3.31 


0.30 


K 



-43- 



Table 2 — Continued 



JD 


RV 


Error 


Observatory 


(-2450000) 


(ms" 1 ) 


(ms- 1 ) 




5024.84232 


2.22 


0.29 


K 


5049.76443 


-2.85 


0.29 


K 


5050.76541 


-1.49 


0.28 


K 


5051.75441 


-1.87 


0.30 


K 


5052.75467 


-0.73 


0.33 


K 


5053.76245 


-1.65 


0.34 


K 



Table 3. Circular Orbit Solutions (Epoch JD 2453369.166) 



Planet 


Period 


M 


e 




K 


Msin i 


a 




(days) 


(degrees) 




(degrees) 


(ms- 1 ) 




(AU) 


61 Vir b 


4.2149±0.0006 


268±13 





n/a 


2.09±0.23 


5.1±0.6 


0.050201±0.000005 


61 Vir c 


38.012±0.036 


210±16 





n/a 


3.58±0.25 


18.2±1.3 


0.2175±0.0001 


61 Vir d 


123.98±0.40 


77±33 





n/a 


3.18±0.29 


24.0±2.2 


0.478±0.001 



Table 4. Keplerian Orbital Solutions (Epoch JD 2453369.166) 



Pla.nct 


Period 


M 


e 


w 


K 


M sini 


a 




(days) 


(degrees) 




(degrees) 


(ms- 1 ) 


(M ffi ) 


(AU) 


61 Vir b 


4.2150±0.0006 


166±53 


0.12±0.11 


105±54 


2.12±0.23 


5.1±0.5 


0.050201±0.000005 


61 Vir c 


38.021±0.034 


177±40 


0.14±0.06 


341±38 


3.62±0.23 


18.2±1.1 


0.2175±0.0001 


61 Vir d 


123.01±0.55 


56±25 


0.35±0.09 


314±20 


3.25±0.39 


22.9±2.6 


0.476±0.001 



